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ABSTRACT 

Aims. Using high-quality, broad-band afterglow data for GRB 091029, we test the validity of the forward-shock model for gamma-ray 
burst afterglows. 

Methods. We used multi-wavelength (NIR to X-ray) follow-up observations obtained with the GROND, BOOTES-3/YA and Stardome 
optical ground-based telescopes, and the UVOT and the XRT onboard the Swift satellite. The resulting data of excellent accuracy allow 
us to construct a multi-wavelength light curve with relative photometric errors as low as 1%, as well as the well-sampled spectral 
energy distribution covering 5 decades in energy. 

Results. The optical/NIR and the X-ray light curves of the afterglow of GRB 091029 are almost totally decoupled. The X-ray light 
curve shows a shallow rise with a peak at ~ 7 ks and a decay slope of a ~ 1.2 afterwards, while the optical/NIR light curve shows 
a much steeper early rise with a peak around 400 s, followed by a shallow decay with temporal index of a ~ 0.6, a bump and a 
steepening of the decay afterwards. The optical/NIR spectral index decreases gradually by over 0.3 before this bump, and then slowly 
increases again, while the X-ray spectral index remains constant throughout the observations. 

Conclusions. To explain the decoupled light curves in the X-ray and optical/NIR domains, a two-component outflow is proposed. 
Several models are tested, including continuous energy injection, components with different electron energy indices and components 
in two different stages of spectral evolution. Only the last model can explain both the decoupled light curves with asynchronous peaks 
and the peculiar SED evolution. However, this model has so many unknown free parameters that we are unable to reliably confirm or 
disprove its validity, making the afterglow of GRB 091029 difficult to explain in the framework of the simplest fireball model. This 
conclusion provides evidence that a scenario beyond the simplistic assumptions is needed to be able to model the growing number of 
well-sampled afterglow light curves. 

Key words, gamma rays: bursts - ISM: jets and outflows - X-rays: individual: GRB 091029 



1. Introduction 

Since the fi rst Gamma-Ray Burst (GRB) was discovered in the 
late 1960's dKlebesadel et alJll973h . the GRB field has evolved 
rapidly, mainly thanks to three generations of dedicated satel- 
lites. The Compton Gamma-Ray Observatory was launched 
in 1991 and with its instrument BATSE showed that GRBs 
are isotropically distributed i n the sky, suggestin g they might 
have a cosmological origin (Meeganetal. 1992). This claim 
was confirmed in 1997 by Beppo-SAX, an Italian-Dutch satel- 
lite that detected and pre cisely localized GRBs at X-ray wave- 
lengths dCosta et al.ll 19971) and enabled ground-based telescopes 

* Present address: American River College, Physics & Astronomy 
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to perform follow-up observations (Ivan Paradiis et alJll997h . in- 
cluding redshift measureme nts that confirmed the cosmologi- 
cal di stances of these events (Met zger et al.ll99"7tlKulkarni et al.l 
1998). Finally, the Swift satellite dGehrels et al.ll2004 . launched 
in 2004, allows for a study of the afterglow emission starting 
very e arly after the GRB is de tected by the Burst Alert Telescope 
(BAT, Barthel mv et al.l l2005). thanks to the rapid slewing capa- 
bility of the spacecraft and a precise localization o f the afterglow 
by on board telescopes sensitive at X-ray (XRT. iBurrows et all 
2005) and ultraviolet/optical (UVOT. lRoming et al.ll2005l) wave- 
lengths. Such precise and early localization allows ground-based 
follow-up telescopes to start observing the afterglow emission 
within tens of seconds of the burst onset. 



2 



R. Filgas et al.: GRB 091029: At the limit of the fireball scenario 



The lea ding model for the afterglow emission is the fire- 
ball model dMeszaros & Reesll997HPiranll999l:lMeszarosl2002L 
IZhang & Meszaros 2004j), where the afterglow arises from the 
synchrotron emission of shock-accelerated electrons in a fire- 
ball interacting with the circum-burst medium. While most of 
the afterglo w light curves prior to Swift were cons istent with 
this model (lHalpernet al.l Il999b IStanek et ail 1 19991) . the more 
recent and detailed light curves of afterglows showed features 
that needed various additions and modifications to the sim- 
plest fir eball model. The ear l y steep decays of th e optical light 
curve s dAkerlof et al.l 11999b iFox et al.l l2003blah iRvkoffet al l 
2004 ) are interpreted as reverse shocks (ISari & Piradll999allbc 
Meszaros & Reesl 119991: iKobavashil bOQOl). Someof the re- 
imps dPar 



brightening s and b umps (Pan aitescu et al 



1998| lGranotetal.1 



2003 - iGuetta et al. l 120071) are attributed to ref reshed shocks 
dRees & Meszaros! 1 19981: lPanaitescuet al.il 19981). the others to 
density variations (Ramirez-Ruiz et al. 2001; 1l)ai&Lull200l 



iets dBercer et al.l 
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2011b). 



However, with the latest generation of GRB instruments capable 
of high sampling in both time and energy domains, the modifi- 
cations made to the standard model still fal l short to explaining 
the observe d afterg lows consistently (e.g., iNardini et al.ll2011t 
lFilgasetal.ll2011at) . 

The Gamma-Ray burst Optic al Near-infrared Detector 
(GROND, iGreiner et aljfeOOSl 120071) has provided high-quality, 
very well-sampled, simultaneous data in seven bands since 
2007, when it was mounted at the 2.2 m MPI/ESO telescope 
at La Silla observatory in Chile. The high-precision data ob- 
tained by GRON D allow for a d e tailed study of afte rglow 
light curv es dGreiner et all l2009bl: INardini et all 1201 ll) . jets 
of GRBs iKrUhler et al l 120091). the dust in their host galax- 
ies dKriihler et al.ll2008HKupcu Yoldas et al.ll2010HGreiner et al.1 
1201 It iKruhler et al.l 120 1 lat ISchadv et al.ll2012l) . their redshifts 
dGreiner et al.l l2009at iKriihler et al.l 1201 lbl) . their associations 
with SNe dOlivares E. et al.l 20121) . and provide tools to test the 



standard fireball scenario and its modifications. 

Here we provide details of the SwiftfXRT, Swift/UVOT, 
GROND, BOOTES-3 and Stardome observations of the after- 
glow of GRB 091029 and discuss the light curves and spectral 
energy distributions (SEDs) in the context of the fireball shock 
model thanks to the very good energy and time-domain cover- 
age of our high-quality data. Throughout the paper, we adopt the 
convention that the flux density of the GRB afterglow can be 
described as F v (t) oc f~ a v~P, where a is the temporal and B the 
spectral index. Unless stated otherwise in the text, all reported 
errors are at 1 <x confidence level. 



2. Observations 

2.1. Swift 

The Swift /BAT was triggered by the long GRB 091029 at 
Tp = 03:53:22 UT and started immediately slewing to the burst 
(Grupe et al. 2009). The mask-weighted light curve shows three 
overlapping peaks, starting at To - 10 s and ending at To + 70 s, 
with peaks at To + 2, +20, and +26 s. The m easured 7*90 (15- 
350 keV) is 39.2 + 5 s dBarfhelmv et al.ll2009h . The BAT prompt 
emission spectrum from To — 1.8 to To + 60.2 s is best fitted 
using a power-law with an exponential cutoff. This fit gives a 
photon index of 1.46 + 0.27 and an £p e ak = 61.4 + 17.5 keV. 
For this model the total flue nce in the 15-150 keV energy range 
is 2.4 + 0.1 x 10" 6 erg/cm 2 dBarthelmv et al.ll2009h . Using stan- 



dard concor dance cosmology (Hp = 71.0 km/s/Mpc, = 0.27, 
Q A = 0.7 3, Komat suetaiTl2009l) . and a redshift of z = 2.752 
(IChornock et al.ll2009l) . the bolometric (IkeV - lOMeV) energy 
release of GRB 091029 is E iso = 8.3xl0 52 erg, with a rest-frame 
Speak of ~230 keV. 

The SwiftfKKY started observ ations of the fie ld of 
GRB 091029 79.3 s after the trigger dGrupe et al.ll2009T) . XRT 
data were obtained from the public Swift archive and reduced in 
the standard manner using the xrtpipeline task from the HEAsoft 
package, with response matrices from the most recent CALDB 
release. The XR T light curve was o btained from the XRT light 
curve repository (Evans et al l2007ll2009l) . 

The Swift/UVOT began settled observa tions of the field o f 
GRB 091029 91 s after the trigger dMarshall & Grupd 12009). 
The afterglow was detected in the White, U, B and V filters. For 
this analysis, we have reduced both image and event mode data. 
Before the count rates were extracted from the event l ists, the as- 
tromet ry was refined following the methodology in Oates et al. 
(2009). The photometry was then extracted from the event lists 
and image files using the FTOOLs uvotevtlc and uvotmaghist, 
respectively, using a source aperture centered on the optical 
position and a background region located in a source-free re- 
gion. We used a 3" source aperture to avoid contamination from 
two neighbouring stars and applied aperture corrections to the 
photo metry in order to be compatible with the UVOT calibra- 
tion dBreeveld et al.ll20TH) . The analysis pipeline used software 
HEADAS 6.10 and UVOT calibration 201 1 1031. In order to be 
compatible with the GROND photometry, UVOT magnitudes 
are provided as AB magnitudes and listed in Table [5] 




Fig. 1. GROND g'-band image of the field of GRB 091029 ob- 
tained 463 s after Tq. The optical afterglow is shown inside the 
Swift XRT error circle with double diameter for better clarity. 
The secondary standard stars are numbered from 1 to 4 and their 
magnitudes reported in Table Q] 



2.2. GROND 

GROND responded to the Swift GRB alert and initiated au- 
tomated observatio ns at 03:57 UT, 4.5 min after the trigger 
(iFilgas etal J 120091) . and imaged the field of GRB 091029 at 
seven later epochs until To + 56 days. A variable point source 
was detected in all bands by the automated GROND pipeline 
(lYoldas et al.l 120081) . The position of the transient was cal- 
culated to be R.A. (J2000) = 04:00:42.62 and Dec (J2000) 
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Table 1. Secondary standards in the GRB 091029 field in the GROND filter bands used for the calibration 
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Star 


R.A., Dec 


g' 


r 




i' 


z' 


number 


[J2000] 


(mag AB ) 


(mag 


ab) 


(mag AB ) 


(mag AB ) 


1 


04:00:47.46, -55:57:35.1 


18.50 + 0.04 


17.03 : 


t0.03 


16.31+0.06 


15.92 + 0.04 


2 


04:00:53.70, -55:57:39.2 


18.63 ± 0.04 


18.14: 


t0.03 


17.98 ± 0.06 


17.85 ±0.05 


3 


04:00:50.68, -55:58:57.3 


19.25 ± 0.04 


18.95 : 


t0.03 


18.92 ±0.07 


18.89 ±0.06 


4 


04:00:48.90, -55:58:43.8 


16.28 + 0.04 


15.75 : 


t0.03 


15.66 + 0.06 


15.55 ±0.04 



Star 


R.A., Dec 


J 




H 


K s 


number 


[J2000] 


(mag Vega ) 


(magvega) 


(mag Vega) 


1 


04:00:37.34, -56:01:20.6 


13.03 ± 


0.03 


12.67 ± 0.03 


12.62 ±0.03 


2 


04:00:39.43, -55:56:02.0 


12.95 ± 


0.03 


12.65 ± 0.03 


12.58 ± 0.03 


3 


04:00:45.75, -55:55:34.7 


13.35 ± 


0.03 


13.07 ± 0.03 


12.96 ± 0.03 


4 


04:00:47.49, -55:57:35.0 


14.66 ± 


0.03 


13.98 ± 0.03 


13.80 ±0.03 



= -55:57:20.0 co mpared to USNO-B reference field stars 
dMonet et al.l l2~003) with an astrometric uncertainty of 0'.'3. 

The optical and NIR image reduct ion and ph otometry were 
performed using standard IRAF tasks (Todv 19931) similar to the 
procedure described in detail in lRriihler et alj (12008). A general 
model for the point-spread function (PSF) of each image was 
constructed using bright field stars and fitted to the afterglow. 
In addition, aperture photometry was carried out, and the results 
were consistent with the reported PSF photometry. All data were 
corrected for a Galactic foregroun d reddening of Er-v = 0.016 
mag in the direction of the burst (Schlege l et aljfl998l) . corre- 
sponding to an extinction of Ay = 0.05 using Ry = 3.1. Optical 
photometric calibration was performed relative to the magni- 
tudes of four secondary standards in the GRB field, shown in 
Fig. [Hand Ta bled] During photometric conditions, an SDSS field 
(Smit h et al.lF2002h at R.A. (J2000) = 03:50:03.25, Dec (J2000) 
= -00:00:37.9 was observed within a few minutes of observa- 
tions of the GRB field. The obtained zeropoints were corrected 
for atmospheric extinction and used to calibrate stars in the 
GRB field. The apparent magnitudes of the afterglow were mea- 
sured with respect to the secondary standards reported in Table 
[TJ The absolute calibration of the JHK S bands was obtained 
with respect to magnitudes of the Two Micron All Sky Survey 
(2MAS S) stars within the GR B field obtained from the 2MASS 
catalog dSkrutskie et al.l l2006) and converted to AB magnitudes. 
All GROND data are listed in Tables|6]and|2] 



2.3. Stardome and BOOTES-3/YA 

The afterglow was observed with the Stardome 0.4 m tele- 
scope, located in Auckland (North Island, New Zealand), us- 
ing a SBIG ST-L-6303E CCD. Images were obtained through 
a OG530 (Schott) filterQthat transmits wavelengths above 5300 
Angstrom. The observations consisted of 300 s exposures that 
were combined in sets of 6 to improve the S/N ratio. Further 
observations were obtained with the Yock-Alle n (YA) robotic 
telesc ope at the BOOTES-3 astronomical station (ICastro-Tiradol 
1201 11) . a fast-slewing 0.6 m Ritchey-Chretien telescope equipped 
with an iXon-889 EMCCD camera located in Blenheim (South 
Island, New Zealand). The observations consisted of a series of 
60 s unfiltered exposures, which were combined in groups to 
improve the S/N ratio. Image reduction was done using standard 
techniques in IRAF. Aperture photometry was performed, using 
PHOT with apertures equivalent to the seeing. For Stardome im- 
ages, a PSF-matching photometry was preferred. All data were 



then cross-calibrated using GROND photometry to obtain con- 
sistent results and are listed in Tables [8] and [9] 



Time since GRB trigger [d] 
0.01 0.1 1 




100 1000 10000 100000 le+06 

Time since GRB trigger [s] 

Fig. 2. Light curve of the X-ray (top panel) and ultraviolet, opti- 
cal and near-infrared (bottom panel) afterglow of GRB 091029. 
Shown data are corrected for Galactic foreground extinction and 
are in AB magnitudes. Upper limits are not shown for better clar- 
ity. Gray regions show the time intervals where XRT data were 
obtained for the broad-band SEDs (Fig.[6]l. 



1 http : //www . optical- filters . co . uk/og538 . html 
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3. Results 

3.1. Afterglow light curve 

The X-ray light curve of the afterglow of GRB 091029 (Fig.© 
shows a very steep decay {a = 3.69 + 0.10) unti l ~ 200 s, con- 
sisten t with being the tail of the GRB emission (Fenimore et al. 
119961) . connecting the prompt phase of the GRB and its after- 
glow. A steep X-ray flare follows after the decay , which de- 
clines rapidly with a = 3.91 + 0.39. Given that the tempo- 
ral decay indices before the flare and after its peak are consis- 
tent within lcr, the possible scenario for this rapid rebrighten- 
ing might involve a r e freshed shock (Rees & Meszaros 1998; 
Panai tescu et all 119981: ISari & Meszarosl l2000t iKumar & Piranl 
2000t IZhang' & Meszarosl 12002b. although th e flare might be 
too rapid for this scenario dHuang et al.1l2006l) . As we are pre- 
dominantly interested in the afterglow phase of the GRB evo- 
lution, we exclude this flare from all our fits. The X-ray light 
curve after ~ 700 s (Fig. |3]l is be st fitted with a broken power- 
law with a smooth break dBeuermann et al.llT999b . The best-fit 
(red. x 1 - 0.89) values of this model are axi = -0.12 + 0.10, 
Wak = 7.4 + 1.8 ks and axi = 1-20 + 0.04. The smoothness s 
iterated to a value of 1, which was set as a lower bound in the 
fit to better constrain the values of the temporal decay indices 
before and after the break. 

Time since GRB trigger [d] 
0.001 0.01 0.1 1 10 
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Fig. 3. The smoothly broken power-law fit to the X-ray light 
curve of the GRB 091029 afterglow. The fitting was applied to 
datapoints after 700 s in order to exclude the early steep transi- 
tion phase and flaring. 

The optical/NIR light curve of the GRB 091029 afterglow 
shows a steep initial rise from the start of the observations un- 
til the peak at around 400 s. The temporal slope of the rise, 
obtained from a fit of UVOT datapoints, is a — -2.90 + 
0.67, consistent with the jet e xpanding in an ISM environment 
(Panaitesc u & Vestrandl20 08). The early peak in the optical/NIR 
light curve is probably not the counterpart of the X-ray flare due 
to the time shift of both peaks. The decay following the initial 
peak has a slope of a = 0.58 + 0.01 until around 5 ks, when it 
starts to flatten. This decay index is obtained from a simultane- 
ous fit to the GROND datapoints between 0.6 - 5 ks in all bands. 

However, fitting this shallow decay phase in each opti- 
cal/NIR band separately shows a steepening of the temporal in- 
dex with increasing wavelength of the GROND filters (see Table 
|2|, suggesting that the afterglow gets bluer. To fit the whole com- 
plex optical/NIR light curve from the beginning to the end of 
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Fig. 4. The two-component fit to the r'-band data obtained by 
GROND, Stardome and BOOTES-3/YA. The parameters of the 
fit for both components (dotted lines) are listed in Table [3] The 
solid line represents the superposition of the two components 
and the best fit to the data. 

the observations, a two-component model is proposed. The first 
component, which dominates the observed optical light curve 
up until the bump at 7q + 5 ks, is composed of three smoothly 
connected power-laws. The second component was needed to 
model the later hump and steep decay, and uses two smoothly 
connected power-laws. The obtained parameters of this fit (red. 
X 1 = 0.92) are listed in Tableland are discussed later. 

Table 2. Light-curve fit parameters for the afterglow of 
GRB 091029 in the time interval of 0.6 - 5 ks. The fitting of the 
NIR bands is affected by the somewhat lower signal-to-noise ra- 
tio of the NIR data as compared to the optical bands. The higher 
value of red. x 1 m the H band is caused by larger residuals. 



Bands 


a 


x 2 /d.o.f. 


g'r'i'z'JHK, 


0.576 ± 0.004 


68/71 


s' 


0.542 ± 0.009 


2.9/ 15 


r 1 


0.574 ± 0.006 


5.1/15 


i' 


0.593 ± 0.010 


3.2/ 15 


z! 


0.622 ±0.018 


6.1/14 


J 


0.601 ± 0.028 


0.8/2 


H 


0.672 ± 0.047 


9.6/2 


K s 


0.815 ±0.075 


2.9/2 



3.2. Afterglow SEDs 

Given that the differences in decay slopes for each GROND fil- 
ter point to a colour evolution, we need to study the SEDs of the 
afterglow. Thanks to the simultaneous multi-band observing ca- 
pabilities of GROND, it is possible to measure the spectral slope 
j8 of the optical/NIR data as a function of time. Fig. [5] shows 
that the optical/NIR spectral index decreases from 0.57 ± 0.04 
to 0.26 + 0.03 between 0.4 and 9 ks, and then slowly increases 
again to a value of 0.49 + 0.12 at around 100 ks. In addition, 
broad-band optical/NIR to X-ray SEDs were constructed at four 
different time intervals within this period, which are indicated 
in the light curve (Fig. Spectra were grouped using the grp- 
pha task and fitted with the GROND data in XSPEC v!2 us- 
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Table 3. Light-curve fit parameters for the whole set of r'-band data obtained by GROND, Stardome and BOOTES-3/YA. 
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Fy(t) 


a, 


tdks] 


Si 




a 2 


t 2 [ks] 


s 2 a 3 


TPL (a) 


-1.95 (fixed) 


0.36 ± 0.02 


2.0 ± 


0.1 


0.84 ± 0.05 


2.80 + 0.15 


2.0 ±0.2 2.83 ±0.20 


DPL <b) 


-0.12 ±0.07 


13.90 ± 0.80 


2.0 ± 


0.1 


1.14 + 0.02 







Smoothly connected triple power-law, describing the narrow jet 
Smoothly connected double power-law, describing the wider jet 



ing x 2 statistics. The combined optical/X-ray SEDs were fitted 
with power-law and broken power-law models and two absorb- 
ing columns: one Galactic foreground with a h ydrogen column 
of N H = 1 ■ 14 x 10 20 cirT 2 ( Kalber la et alj2005l) and another one 
that is local to the GRB host galaxy at z = 2.75. Only the latter 
was allowed to vary in the fits. To investigate the dust redden- 
ing in the GRB environment, the zdust model was used, which 
contains Large and Small Magellanic Cl oud s (LMC , SMC) and 
Milky Way (MW) extinction laws from|Pel d!992l) . Fits of op- 
tical/NIR data alone as well as the broad-band fits resulted in 
a host dust extinction that was consistent with zero, therefore in 
all the models we assumed no host dust extinction for simplicity. 
With photometric data alone it was n ot possible to constrain the 
presence of Lyman-alpha absorption dLamb & Reichartl 2000) in 
the g' band from neutral hydrogen within the host galaxy. The 
g'-band data were therefore removed from the SED fits. 



Wavelength [nm] 
100 10 




Time since trigger [s] 

Fig. 5. The optical/NIR spectral slope as a function of time. 

Given that the broad-band SEDs proved to be inconsistent 
with a simple power-law model (red. x 2 - 16.5), we used mod- 
els that include a break between the X-ray and optical/NIR data. 
We fitted all four epochs of broad-band SEDs simultaneously 
with a sharp broken power-law model, where the host-intrinsic 
absorbing column density and the X-ray spectral index are tied 
between each SED but left free to vary (Fig. [6]), due to the fact 
that the X-ray hardness ratio does not change during the after- 
glow. The low-energy spectral indices and energy of the break 
were left untied between SEDs and free to vary. The best fit (red. 
X 2 = 0.94) gives values of the host-equivalent neutral hydrogen 
density Nh,x - (3.0 + 1.3) x 10 21 crrT 2 and a high-energy spec- 
tral index fix = 1-08*q[^. The best-fit values of the low-energy 
spectral indices and break energies are listed in Table [4] This 
fit shows that the break evolves in time to shorter wavelengths 




Fig. 6. Broad-band optical/NIR to X-ray SEDs fitted with a bro- 
ken power-law. Mid-times of each SED are listed in Tableland 
regions where the X-ray data were taken for each SED are shown 
in Fig. [2] 

between SEDs I and II, and then it evolves the opposite way to- 
wards longer wavelengths between SEDs II and IV. It also shows 
that below and above the cooling break A/3 + 0.5 (with quite 
high an d variable significance), which is a value demanded by 
theory dSari et alj|1998l) . This is another indication that we are 
not seeing a simple single radiating electron population. 

Table 4. Best-fit parameters resulting from the broken power- 
law fit to the broad-band SEDs. The host-equivalent neutral hy- 
drogen density N H = (3.0 + 1.3) x 10 21 ctrT 2 . 



SED 
number 


Midtime [s] 
of SED 


Low energy 
spectral index 


Cooling 
break [eV] 


High energy 
spectral index 


I 

11 
111 

IV 


883 
7161 
18056 
104026 


46 +0 06 
32 +0 05 
34 +0 06 

"""-0.06 

42 +0 23 


26 4 +15 - 2 
47 7+20.8 

40 4 +1M 

31 6 +1155 
:l -no 


1.083* 



3.3. Closure relations 

Using values obtained from the diffe rent 
the clo sure relations ([Granot & Saril 

2004 



2001; Zhang & Meszaros 



Pan aitescu & Kumar] |2000t) between 
indices. Assuming the break in the broad-band SEDs to be 



fits we can test 
2002; iDai & C heng 
acusin etlril [2009; 
and spectral 



temporal 
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the cooling break v c , we see that the X-ray data are above this 
frequency and the optical/NIR data below it during the afterglow 
observations. The fit-derived X-ray spectral index fix = 1.08^S"£5 
results in a power-law index of the electron energy distribution 
p = 2.17 ± 0.11. This spectral index and the late X-ray temporal 
slope of axi - 1-20 ± 0-04 a re within Icr consistent with the 
equation (Racusinetal. 2009) for p > 2 and a constant decay 
in the v > v c regime, where the jet is interacting with either a 
homogeneous interstellar medium (ISM) or a wind and is in the 
slow or fast cooling phase. 

The late (f > 20 ks) optical/NIR single-component model 
decay index of a = 1.14 + 0.02 and the late spectral index of 
f3 = 0.49 + 0.12 are within l<x consistent with the equation 
for a normal decay in the v,„ < v < v c regime, where the jet 
is interacting with a wind medium and is in the slow cooling 
phase. However, the steep initial onset of the optical/NIR after- 
glow seems to exclude this scenario, as it is well consistent with 
the expected temporal index of a ~ -3 for the jet interacting 
with the ISM and v c > v opt , wh ereas the fastest possible ris e 
for a wind medium is a ~ -0.5 dPanaitescu & Vestrandl f2008). 
Moreover, during the early shallow decay in the optical/NIR 
bands, the spectral index is evolving with time and thus can- 
not be tested with the simple closure relations. In the case of the 
two-component scenario (Fig. |4|, the resulting spectral index /3 
is the superposition of the spectral indices of the narrow and the 
wide jet and evolves with time as the ratio between these two jets 
changes. Without knowing the spectral indices of each compo- 
nent, the closure relations for such a scenario cannot be tested. 

4. Discussion 

The almost total decoupling of the optical/NIR and the X- 
ray light curves of the afterglow of GRB 091029 suggests 
a double outflow origin (Huang et ail 120041: iPeng et alJ 1 2005: 
Shet hetal] 120031: Uin et all 120071: IWu et all 120051: i Resmi et alj 
2005t iRacusin et alJ 1200a iFilgas et al.1 1201 lbh . This is sup- 
ported by our finding that the X-ray spectral hardness does 
not evolve synchronously with the optical spectral hardening at 
0.3 - 10 ks. We discuss three possible scenarios, all based on a 
two-component jet, to explain the peculiar behavior of the light 
curve and the spectrum of this afterglow. 

4. 1 . Continuous energy injection 

In this scenario, the X-ray light curve after 700 s is produced 
by the wider, X-ray-dominating outflow, which has a decelera- 
tion time of a few ks. The pre-deceleration phase of the wide 
jet would cause the early shallow rise of the X-ray afterglow 
with the peak around 7.4 ks, after which the wide jet would 
turn into normal deceleration producing the axi ~ 1.2 decay. 
The same principle would apply to the early optical light curve, 
where the early rise and peak at ~ 400 s would be a result of 
a pre-deceleration phase of the narrower, optically dominating 
outflow. From the peak time of this narrower outflow, we can es- 
timate Jh^initialLorentz factor in case of the ISM to be T n ~ 300 
using Molinari et a il (120071) or T„ ~ 120 using Ghirlanda et alJ 
d2012l) . in both cases substituting the measured quantities and 
normalizing to the typical values n — 1 cirT 3 and rj - 0.2 
dBloom et alJl2003h . The shallow decay of the optical/NIR light 
curve until 10 ks is then a result of so me form of a continuous 
energy injection by the central engine dRees & Meszaros 1998; 
Dai & Lufll998t iPanaitescu et al.lll998l: ISari & Meszarosp OOO: 
Zhang & Meszarosll2001l: IZhang et al.ll2006t iRossi et aklboi lh . 
When this energy injection in the narrow jet ends at ~ 20 ks, the 



temporal slope steepens to a ~ 1.1, characteristic for a normal 
decay. However, this scenario cannot explain the time evolution 
of the optical/NIR spectral index during the injection time inter- 
val. The hardening of the optical spectrum would require that the 
electron index p changes in the outflow with time and after the 
injection ends, it changes back close to its original value (Fig. 
01, making this scenario somewhat contrived. 

Another issue with this scenario is the fact that in the stan- 
dard interpretations, the narrow jet dominates the X-ray after- 
glow as it is more energetic, and the wider jet dominates the 
optical domain, especially at later times. The X-ray light curve 
is therefore expected to peak before the optical one. In such a 
standard interpretation, the reverse order of peaks in the light 
curves of the afterglow of GRB 091029 could be the result of an 
off-axis geometry of the jet as seen by the observer. If the ob- 
server was inside the cone of the wide, optically dominating jet, 
but outside of the beaming cone of the narrow, X-ray dominating 
jet, the optical light curve could start with a steep rise and early 
peak due to the pre-deceleration phase of the wide jet, and the 
X-ray light curve could follow with a later peak as the relativis- 
tically beamed emission cone of the narrow jet widens and grad- 
ually enters the sightline of the observer. However, given that 
neither the early shallow decay of the optical/NIR light curve 
nor the colour evolution in the same domain are correlated with 
the temporal evolution of the narrow jet, this off-axis scenario 
has the same problems with the unfeasible temporal evolution of 
the electron energy distribution as the previous one. 

4.2. Two outflows with different p parameters 

The second scenario uses a different two-component model, 
without energy injection, where the optical/NIR light curve is a 
superposition of two components as described in section 3.1 and 
shown in Fig. [4] The first component would represent a narrow, 
ultra-relativistic jet, with deceleration time of ~ 400 s, a normal 
decay phase afterwards, and a jet break at ~ 2.8 ks, followed 
by a steep post-jet-break decay. The second component repre- 
sents a wider and mildly relativistic jet with the deceleration 
time of ~ 14 ks and a normal decay. While the residuals in Fig. [4] 
might hint at a possible jet break of the wider component around 
100 ks, the data are within 3cr of the straight power-law and there 
is no visible break in the X-ray data at that time. Therefore we 
cannot make a reliable statement about the presence of such a 
jet break. The narrow jet would be dominant in the early part of 
the light curve and the wide one would be responsible for the 
late hump and dominate the optical/NIR light curve afterwards. 
The shallow decay between 0.4 - 5 ks would be a result of the 
superposition of fluxes from both outflows. Since the emission 
from the narrow component in the normal decay phase is de- 
creasing and the the emission from the wide component in the 
pre-deceleration phase is increasing with time, the ratio of the 
fluxes of the two outflows in this time interval would vary. This 
model can explain the flattening of the spectral index during this 
period, assuming that each outflow with different Lorentz factors 
has a different electron energy distribution index p. As the ratio 
between the narrow jet with a higher/? value (p ~ 2.1, estimated 
from the highest value of the optical/NIR spectral slope) and the 
wide jet with a lower p value (p ~ 1 .5, estimated from the lowest 
value of the optical/NIR spectral slope) changes, the spectral in- 
dex of the co-added flux is gradually evolving as well. The spec- 
tral index from the start of the GROND data is governed by the 
narrow outflow with the soft spectrum but gradually decreases 
as the flux from the wide component with the harder spectrum 
gets dominant (Fig. @). However, this model does not explain 
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the turnover of the spectral evolution at ~ 10 ks, where only the 
wide component with a constant p is dominant. The model also 
does not fit the shape of the X-ray light curve, which should be 
dominated the whole time by the narrow jet due to the X-ray 
spectral index fix being constant and consistent with p ~ 2.1 of 
the narrow jet. 

4.3. Two outflows in different stages of the spectral evolution 

The third scenario uses a slightly modified two-component jet 
setup described in section 4.2 (Fig. 0), in which now both out- 
flows have the same value of p ~ 2. 1 and are at different stages 
of the synchrotron spectral evolution dSari et al.H l998). Both the 
X-ray and optical/NIR light curves can be modeled as superpo- 
sitions of the two components. The flattening of the SEDs II and 
III in the optical/NIR region (Fig. |6]l would then be a result of 
the wide jet having both the cooling break v c and the injection 
frequency v m between the X-ra y and optical/NIR wavelengths 
dNicuesa Guelbenzu et al.ll2012l) . while the narrow jet has only 
the cooling frequency v c between X-ray and optical/NIR bands. 
In that case, the spectral slope in the optical/NIR bands of the 
narrow jet woul d befe-0.5 = .58±0.06 and that of the wide jet 
would be -1/3 dSari et al.lll998h . As the ratio between these two 
outflows changes, it would explain the spectral hardening in the 
optical/NIR bands, while being consistent with the X-ray spec- 
tral slope staying constant thanks to equal p values in both out- 
flows. The turnover in the spectral evolution at ~ 10 ks can be ex- 
plained by the passage of the frequency v m through the GROND 
filters, after which the optical/NIR spectral index would be con- 
sistent with the spectral phase of the narrow jet. Given that the 
softening of the optical/NIR spectrum after ~ 10 ks is slow and 
gradual, the spectral break at the frequency v„ must be very 
smooth dGranot & Sarill2002tlFilgas et alj|201 lab . 

This model can also explain different break times in the X- 
ray and optical/NIR light curves, assuming that the deceleration 
time of the wide jet is defined by the X-ray light-curve peak but 
the movement of the break v m is counter-balancing the expected 
decay steepening in the optical/NIR bands until it passes through 
them and causes the late break in the light curve at Tq = 10 ks. 
Of course, it is questionable whether the movement of the break 
v m could counter-balance the turnover visible in the X-ray light 
curve so perfectly that it would in fact completely negate it in 
the optical/NIR domain and produce such a straight power-law 
as we see in the optical/NIR light curve between 0.4 - 5 ks. 
This scenario is difficult to confirm or disprove, though, by fit- 
ting the light curve and SEDs alone because this model has a 
large number of free parameters. Therefore it is very difficult 
to fit the broad-band SEDs with a model consisting of a double 
power-law and a triple power-law component and constrain all 
the spectral slopes and synchrotron break frequencies reliably. 
The fitting of light curves is confronted with similar difficulties. 
The optical/NIR light curve might need an even more complex 
model than the one presented in Fig. [4] in order to describe the 
effects of the moving frequency v m . However, as our data are the 
result of the superposition of two components, fitting of the data 
does not provide us with reliable constraints on the parameters 
of each component. 

5. Conclusions 

The growing number of well-sampled data sets from the latest 
generation of instruments like the Swift satellite and the GROND 
imager show that the radiative mechanism responsible for the 
optical to X-ray GRB afterglow emission is not as simple and 



well understood as previously believed. The simplest fireball 
model has an increasingly difficult time to explain the complex 
light curves of some GRB afterglows. In many cases, the opti- 
cal and X-ray emission are seemingly decoupled, thus providing 
an indication that they are produced by different mechanisms. 
The afterglow of GRB 091029 is an extreme case, where the 
optical/NIR and the X-ray light curves are almost totally decou- 
pled, as if they belonged to two different GRBs. Moreover, the 
GROND SEDs show a strong colour evolution with the opti- 
cal/NIR spectral index decreasing from 0.57 to 0.26 between 0.4 
and 9 ks, and then increasing again to a value of ~ 0.49 at around 
100 ks, while the X-ray spectral index remains constant through- 
out the observations. This observational evidence leads us to the 
conclusion that the emission in both energy bands needs to be 
produced by two different outflows. 

We discuss several possible scenarios to explain this peculiar 
afterglow. The first one includes a continuous energy injection to 
explain the shallow initial decay of the optical/NIR light curve. 
However, this model is not able to explain the spectral evolu- 
tion during the injection period, given that the theory assumes 
the electron energy distribution index p of the outflow is con- 
stant. To solve this, the second scenario uses two components 
with different p values. As the ratio between these two outflows 
changes, the resulting spectral index changes as well. This model 
is, however, not able to explain either the turnover in the opti- 
cal/NIR spectral-slope evolution, nor the different times of the 
breaks in the X-ray and optical/NIR light curves. The third sce- 
nario offers a solution by putting the two outflows with similar 
p values into two different stages of the spectral evolution. The 
narrow jet, dominating the optical/NIR wavelengths before the 
hump, has a cooling break between the optical and the X-ray 
bands, while the wide jet, responsible for the late optical/NIR 
light curve, has both v c and v m frequencies between the optical 
and the X-ray bands. During the hump, the injection frequency 
v m passes through the GROND filters and the light curve be- 
comes similar to the X-ray one. 

Even though the last model can in principle explain the irreg- 
ularities in the afterglow of GRB 091029, its complexity does 
not allow us to test it reliably, despite the large high-quality 
dataset presented in this work. Ironically, the data quality of 
the presented afterglow light curve is so good it rules out any 
simple model for the temporal-spectral evolution, but is not 
good enough to really constrain the more complex, constructed 
models. Rather than using the forward-shock scenario, alterna- 
tive models might be needed to explain the multi-wavelength 
data of the afterglow of GRB 091029. For example a reverse- 
forward shock emission, where th e optical afterglow is domi- 
nated by a long-li ved reverse shock dUh m & Beloborodov 120071: 
iGenet et ani2007l) and the X-ray afterglow is from the forward 
shock. This model can decouple the two light curves almost 
completely, requires mass injection into the blast-wave and has 
many free parameters. Another alternative can be the cannon- 
ball model (Dado & Darll20loh . where the observed afterglow 
emission is described as the sum of thermal bremsstrahlung and 
synchrotron emission produced by one or more cannonballs de- 
celerating in the circ umburst medium. The last we mention is t he 
late prompt model (|Ghisellini et ai]|2007l: iKumar et alJl2008albt 
Ghis ellini et al.ll2009t iNa rdini et alJ l2010h . which was proposed 
to explain the different temporal and spectral evolution in the 
optical and X-ray bands of certain long GRBs which show a 
late-time chromatic flattening. The authors interpret the complex 
broadband evolution as due to the sum of two separate processes: 
the standard forward shock and the emission produced by a late- 
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time activity of the central engine (i.e., the so called late prompt 
emission). 
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Table 5. UVOT photometric data 



T m id - To [ks] 


Exposure [s] 


Brightness^' magAB 


Filter 


0.1156 


50 


22.18 + 0.71 


White 


0.1656 


50 


21.39 + 0.35 


White 


0.2155 


50 


20.30 + 0.16 


White 


0.5932 


20 


19.45 + 0.16 


White 


0.7674 


20 


19.71+0.19 


White 


0.9332 


150 


19.68 + 0.08 


White 


1.2584 


193 


19.81+0.14 


White 


1.6049 


193 


19.81+0.14 


White 


2.0387 


366 


20.23 + 0.14 


White 


2.3850 


20 


19.94 + 0.21 


White 


6.6930 


200 


20.60 + 0.11 


White 


8.1318 


200 


20.77 + 0.12 


White 


50.5180 


6487 


22.20 + 0.13 


White 


65.1696 


12382 


22.32 + 0.13 


White 


122.7112 


23679 


23.00 + 0.24 


White 


146.0998 


12138 


> 22.94 


White 


203.9169 


35415 


> 23.06 


White 


311.1750 


40962 


> 23.11 


White 


404.1329 


75705 


> 23.4 


White 


496.6075 


98012 


> 23.91 


White 


554.1967 


6354 


> 22.91 


White 


0.3277 


50 


19.79 + 0.25 


U 


0.3777 


50 


19.36 + 0.20 


U 


0.4277 


50 


19.09 + 0.17 


U 


0.4777 


50 


19.64 + 0.23 


u 


0.5276 


50 


19.29 + 0.19 


u 


0.7181 


20 


20.23 + 0.52 


u 


1.2085 


193 


20.35 + 0.40 


u 


1.5553 


193 


20.29 + 0.38 


u 


1.9896 


366 


20.00 + 0.26 


u 


2.3351 


20 


19.79 + 0.40 


u 


6.2829 


200 


20.96 + 0.26 


u 


7.7213 


200 


20.91+0.25 


u 


18.9286 


908 


21.42 + 0.17 


u 


25.3918 


444 


21.77 + 0.32 


u 


36.9765 


414 


> 21.93 


u 


49.8063 


6486 


> 22.07 


u 


64.5991 


12538 


22.25 + 0.24 


u 


75.6289 


337 


> 21.39 


u 


122.3652 


23602 


> 22.96 


u 


145.6560 


12214 


> 22.10 


u 


203.5215 


35436 


> 22.73 


u 


310.5329 


40872 


> 23.20 


u 


403.6852 


75885 


> 22.98 


u 


473.2423 


52055 


> 22.84 


u 


635.9086 


75370 


> 23.60 


u 


985.0417 


47426 


> 23.42 


u 


1340.8648 


52778 


> 23.74 


u 


1669.9707 


88045 


> 23.81 


u 


2034.5151 


41876 


> 23.66 


u 


0.5689 


20 


18.84 + 0.32 


B 


0.7431 


20 


18.50 + 0.26 


B 


1.2336 


194 


18.88 + 0.23 


B 


1.5803 


193 


18.72 + 0.20 


B 


2.0143 


366 


18.92 + 0.19 


B 


2.3598 


20 


19.51+0.51 


B 


6.4882 


200 


19.63 + 0.17 


B 


7.9271 


200 


19.63 + 0.17 


B 


19.6085 


441 


20.15 + 0.16 


B 


310.8542 


40917 


> 21.78 


B 


403.9090 


75794 


> 22.88 


B 
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Table 5. continued. 



R. Filgas et al.: GRB 091029: At the limit of the fireball scenario 



T m id - T [ks] 


Exposure [s] 


Brightness^' mag A B 


Filter 


473.4810 


52145 


> 23.00 


B 


0.0745 


10 


> 18.19 


V 


0.6446 


20 


17.70 ± 0.29 


V 


0.8170 


20 


> 17.85 


V 


1.1348 


193 


18.33 + 0.29 


V 


1.4813 


192 


18.58 + 0.34 


V 


1.9159 


366 


> 18.59 


V 


2.3483 


193 


18.82 + 0.40 


V 


7.1045 


200 


19.30 + 0.25 


V 


13.1297 


907 


19.25 + 0.11 


V 


30.4949 


907 


20.08 ± 0.20 


V 


51.2643 


6556 


20.72 + 0.26 


V 


65. /636 


122/0 


20. / / + 0.25 


T / 

V 


123.0766 


23794 


> 21.18 


V 


146.5634 


12101 


> 21.11 


V 


204.3240 


35420 


> 21.77 


V 


311.4738 


40960 


> 21.01 


V 


407.2422 


81383 


> 20.92 


V 


476.8197 


46508 


> 20.98 


V 



Corrected for Galactic foreground reddening. 
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Tmid - T Q [ks] 


Exposure [s] 


8' 


Brightness'"' mag AB 
r' i' 


z' 


0.3067 


66 


18.19 + 0.04 


17.84 ±0.02 


17.78 ±0.03 


17.59 ±0.08 


0.4103 


66 


18.00 ± 0.04 


17.69 ± 0.02 


17.60 ± 0.02 


17.38 ± 0.04 


0.5154 


66 


18.09 ± 0.03 


17.76 ± 0.02 


17.66 ± 0.03 


17.49 ± 0.05 


0.7290 


66 


18.22 ± 0.03 


17.92 ± 0.02 


17.87 ± 0.04 


17.64 ± 0.05 


0.8257 


66 


18.31 ±0.03 


17.98 ± 0.02 


17.94 ± 0.03 


17.76 ± 0.06 


0.9237 


66 


18.39 ± 0.03 


18.07 ± 0.02 


18.02 ±0.03 


17.84 ±0.05 


1.0238 


66 


18.47 ± 0.04 


18.14 ±0.02 


18.09 ±0.03 


17.92 ± 0.07 


1.1549 


115 


18.50 ± 0.03 


18.21 ±0.02 


18.16 ±0.04 


18.02 ± 0.05 


1.3407 


115 


18.58 ± 0.04 


18.31 ±0.01 


18.26 ±0.03 


18.12 ±0.05 


1.5249 


115 


18.67 ± 0.03 


18.39 ± 0.02 


18.36 ±0.02 


18.19 ±0.05 


1.7169 


115 


18.73 ± 0.02 


18.46 ±0.01 


18.44 ±0.03 


- 


1.9195 


115 


18.80 ± 0.03 


18.53 ±0.02 


18.50 ±0.03 


18.32 ±0.05 


2.1132 


115 


18.83 ± 0.03 


18.58 ±0.02 


18.56 ±0.02 


18.44 ± 0.04 


2.3095 


115 


18.90 ± 0.03 


18.64 ±0.01 


18.63 ±0.03 


18.43 ± 0.05 


2.5040 


115 


18.96 ± 0.03 


18.70 ± 0.01 


18.66 ±0.02 


18.51 ±0.05 


2.8368 


375 


19.04 ± 0.02 


18.77 ±0.01 


18.76 ±0.02 


18.58 ±0.03 


3.2896 


375 


19.12 ±0.02 


18.88 ±0.01 


18.85 ±0.01 


18.73 ±0.03 


3.7418 


375 


19.20 ± 0.02 


18.95 ± 0.02 


18.93 ±0.02 


18.81 ±0.03 


4.1939 


375 


19.27 ± 0.02 


19.02 ± 0.01 


19.00 ± 0.02 


18.86 ±0.04 


4.6568 


375 


19.32 ± 0.02 


19.07 ±0.01 


19.05 ± 0.02 


18.89 ± 0.03 


5.1021 


375 


19.35 ± 0.02 


19.11 ±0.01 


19.09 ± 0.02 


18.98 ± 0.04 


5.5542 


375 


19.38 ±0.01 


19.14 ±0.01 


19.11 ±0.02 


18.98 ±0.03 


6.0063 


375 


19.39 ± 0.02 


19.15 ±0.01 


19.14 ±0.02 


19.02 ± 0.04 


6.4708 


375 


19.40 ± 0.02 


19.17 ±0.01 


19.15 ±0.02 


19.01 ±0.03 


6.9143 


375 


19.41 ± 0.02 


19.18 ±0.01 


19.17 ±0.02 


19.04 ± 0.03 


7.3624 


375 


19.42 ± 0.01 


19.18 ±0.01 


19.17 ±0.02 


19.06 ± 0.03 


7.8060 


375 


19.45 ± 0.01 


19.21 ±0.01 


19.18 ±0.02 


19.05 ± 0.03 


8.2685 


375 


19.46 ± 0.02 


19.21 ±0.01 


19.20 ± 0.02 


19.07 ± 0.04 


8.7202 


375 


19.48 ± 0.02 


19.24 ±0.01 


19.23 ± 0.02 


19.09 ± 0.03 


9.1693 


375 


19.50 ± 0.02 


19.27 ±0.01 


19.26 ± 0.02 


19.10 ±0.04 


9.6214 


375 


19.53 ±0.01 


19.29 ± 0.01 


19.29 ± 0.02 


19.14 ±0.03 


14.4446 


375 


19.70 ±0.01 


19.47 ±0.01 


19.43 ± 0.02 


19.30 ± 0.04 


15.3777 


375 


19.75 ± 0.01 


19.50 ±0.01 


19.48 ± 0.02 


19.32 ± 0.05 


15.8274 


375 


19.78 ± 0.01 


19.53 ±0.01 


19.51 ±0.02 


19.38 ± 0.03 


16.2785 


375 


19.80 ± 0.01 


19.54 ± 0.01 


19.52 ± 0.02 


19.39 ± 0.04 


16.7300 


375 


19.81 ±0.01 


19.57 ±0.01 


19.50 ± 0.02 


19.36 ± 0.04 


17.3570 


375 


19.84 ± 0.02 


19.58 ±0.01 


19.55 ±0.02 


19.41 ±0.04 


17.8057 


375 


19.87 ± 0.02 


19.60 ±0.01 


19.58 ±0.02 


19.42 ± 0.04 


18.2565 


375 


19.95 ± 0.03 


19.63 ± 0.01 


19.55 ± 0.04 


19.46 ±0.05 


104.0014 


686 


21.82 ±0.05 


21.51 ±0.03 


21.44 ±0.06 


21.25 ±0.09 


170.5107 


1714 


22.47 ± 0.09 


22.14 ±0.06 


22.19 ±0.09 


21.95 ±0.11 


188.6865 


1724 


22.62 ± 0.03 


22.34 ± 0.03 


22.43 ± 0.06 


22.12 ±0.09 


344.9195 


3521 


23.28 ± 0.20 


23.21 ±0.15 


23.12 ±0.15 


23.29 ±0.17 


438.7789 


3520 


> 23.83 


23.80 ± 0.23 


23.64 ± 0.28 


> 23.48 


872.9788 


7136 


> 24.19 


> 23.92 


> 23.87 


> 23.85 


1478.9968 


7096 


> 24.04 


> 24.24 


> 23.88 


> 24.10 


4832.0789 


7182 


> 24.81 


> 25.50 


> 24.02 


> 24.37 



Corrected for Galactic foreground reddening. 
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Table 7. GROND JHK S photometric data 



7mid - To [ks] Exposure [s] Brightness' 3 ' mag' 







J 


H 


K s 


0.3135 


82 


17.33 +0.08 


17.26 ± 0.06 


17.17 ± 0.11 


0.4171 


82 


17.22 ±0.08 


17.13 ± 0.06 


16.91 ± 0.09 


0.5221 


82 


17.34 ± 0.07 


17.23 ± 0.06 


17.13 ±0.11 


0.8835 


377 


17.72 + 0.03 


17.61 ± 0.05 


17.28 ± 0.07 


1.4625 


729 


18.09 ±0.03 


17.94 ± 0.04 


17.59 ± 0.08 


2.2379 


754 


18.35 ±0.03 


18.45 ± 0.06 


17.98 ±0.10 


3.5383 


1780 


18.63 ±0.03 


18.57 ± 0.05 


18.55 ±0.11 


5.3549 


1772 


18.91 ±0.04 


18.80 ±0.05 


18.68 ±0.13 


7.1614 


1758 


18.98 ±0.04 


18.86 ± 0.06 


18.68 ±0.11 


8.9684 


1775 


19.01 ± 0.04 


19.12 ± 0.06 


18.91 ±0.14 


16.0777 


1777 


19.31 ±0.05 


19.24 ± 0.08 


18.98 ±0.17 


18.0561 


1775 


19.33 ±0.05 


19.30 ± 0.07 


19.11 ±0.19 


104.0267 


739 


21.33 ±0.40 


> 20.16 


> 19.53 


105.2524 


1644 


21.27 ±0.29 


> 20.48 


> 19.70 


170.5336 


1762 


21.79 ±0.38 


> 20.83 


> 19.98 


188.7100 


1773 


21.90 ±0.32 


> 21.12 


> 20.23 


344.9423 


3569 


> 21.77 


> 21.40 


> 20.66 


438.8024 


3569 


> 21.68 


> 21.10 


> 20.80 


873.0025 


7184 


> 22.03 


> 21.43 


> 21.14 


1479.0224 


7142 


> 22.07 


> 21.51 


> 20.98 


4832.1018 


7230 


> 22.12 


> 21.54 





Corrected for Galactic foreground reddening. Converted to AB magnitudes for consistency with Table[6] 

For the SED fitting, the additional error of the absolute calibration of 0.05 (J and H) and 0.07 (K s ) mag was added. 



Table 8. BOOTES-3/YA photometric data 



Tmid - T [ks] 


Exposure [s] 


Brightness* 


a) magAB 


24.8866 


1440 


20.02 ± 


0.11 


27.4346 


1980 


20.09 ± 


0.15 


30.1346 


1740 


20.02 ± 


0.12 


31.9248 


1800 


20.10 ± 


0.10 


36.4357 


4320 


20.22 ± 


0.10 


42.6773 


5220 


20.46 ± 


0.12 



(a) Corrected for Galactic foreground reddening. 
Table 9. Stardome photometric data 



?mid - T [ks] 


Exposure [s] 


Brightness* 


a) mag AB 


14.2620 


1800 


19.42 ± 


0.20 


16.1568 


1800 


19.56 ± 


0.14 


18.0697 


1800 


19.48 ± 


0.14 


19.9627 


1800 


19.51 ± 


0.16 


21.8566 


1800 


19.82 ± 


0.16 


23.4610 


1800 


19.70 ± 


0.14 


25.6228 


1800 


20.10 ± 


0.14 


30.0153 


1800 


20.13 ± 


0.16 


31.9317 


1800 


20.11 ± 


0.14 


34.4200 


1800 


20.31 ± 


0.12 


37.3110 


1800 


20.26 ± 


0.14 


39.2377 


1800 


20.35 ± 


0.16 


41.1670 


1800 


20.69 ± 


0.15 


43.0920 


1800 


20.64 ± 


0.13 



Corrected for Galactic foreground reddening. 



